Improved Lattice Actions with Chemical Potential 
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We give a prescription how to include a chemical potential fi into a general lattice 
action. This inclusion does not cause any lattice artifacts. Hence its application to an 
improved - or even perfect - action at /i = yields an improved resp. perfect action at 
arbitrary /i. For short-ranged improved actions, a good scaling behavior holds over a wide 
region, and the upper bound for the baryon density - which is known for the standard 
lattice actions - can be exceeded. 

1. Introduction 

At chemical potential /i 7^ 0, the action of QCD is complex (see subtitle of this work- 
shop). With respect to lattice simulations, this means that the standard Monte Carlo 
techniques fail. In principle, it is possible to simulate at /x = and include the baryon 
density hb in the measured observables by a suitable re-weighting. However, in practice 
this is extremely tedious: the ratio of an observable at /i 7^ divided by the observable 
at /i = is exponentially suppressed by the physical volume. Hence the absolute values 
of the statistical contributions tend to be many orders of magnitude larger than the ratio 
of interest, so that tremendous statistics are required ("sign problem") [] 

In addition to the statistical error, lattice simulations are also plagued by systematic 
errors. The worst source of them are artifacts due to the finite lattice spacing a > 0. 
Improved lattice actions are designed to suppress these artifacts, i.e. the continuum scaling 
should persist to a good approximation down to a rather short correlation length in lattice 
units, ^/a, or up to rather large values of /xa (below we use lattice units, a = 1). 

As a particularly bad manifestation of lattice artifacts, there is an upper bound for ub, 
resp. for the fermion density Uf. We will show that improved actions can weaken this 
unphysical saturation effect. 

In simulations, chiral symmetry appears to be restored already at half the pion mass 
("onset problem" 0J§]). Usually quenching was blamed for that, but attempts to go 
beyond the quenched approximation have not really helped in this respect. The Glasgow 
group now assumes that this is an effect of low statistics, but it is also conceivable that 
lattice artifacts contribute to this problem. If this is true, then the use of dynamical 
fermions together with an improved action should help to obtain a more continuum-like 
chiral behavior. 

The extreme case of improved actions are perfect actions: they eliminate all lattice 
artifacts. Unfortunately they tend to involve an infinite number of couplings, hence we 

^In principle, the complex Langevin algorithm is an alternative [P, but it has not lead to QCD results, 
again due to the limited statistics. 
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can only apply short-ranged approximations. A perfect action is constructed by block 
variable renormalization group transformations (RGTs) of the type 

g-5'[f',*',c/'] ^ J jj^jj^jju e-s[^,^,u] exp{-r[^', ^, ^, U', U]} . (1) 

We start from a fine lattice with action S[^, U], and introduce a new lattice with a n 
times coarser lattice spacing, where the action is given by S'[^', \E'', U'] according to eq. 
(|1]). The transformation term T relates the original lattice variables to the new ones by 
averaging in some way over the fine lattice sites close to a given coarse site. It has to 
provide the invariance of the partition functions, Z' = Z, and of the expectation values, 
hence of the physical contents of the theory. The choice for T characterizes the RGT. After 
the RGT, the correlation length in lattice units is reduced by a factor n, C^'/a' = 1/n-^/a. 

Assume a mass m and an inverse temperature (3. Now we start from a very fine lattice 
at mass m/nN and inverse temperature (5nN , and we perform RGTs of block factor 
n. This leads to the desired parameters m, (3. In the limit nN oo, this method yields 
a perfect action at the desired finite parameters m, (3. 

2. Inclusion of the chemical potential /i 

In the continuum, the chemical potential can be included by the substitutions ip{—p) 
i'i—p, —p4 + ifi), i'ip) — ^ i'{p,P'i — ijj) in the quark fields, or by replacing p4 ^ p4 + iii 
in the Dirac operator. It is not trivial how to apply this substitution to the naive lattice 
Dirac operator i'-^^smpy + m. An early guess, sinp4 — sinj94 + i/i does not have the 
correct continuum limit. Instead one should use sin(p4 + i^) Also for the somewhat 
more complicated Wilson-Dirac operator, i'-^ysmpu + m+ {r /2)Y,uPii ipv = 2 sin(pi,/2)), 
the substitution p4 ^ p4 + z/i works. 

If we want to construct a perfect action, we can use a standard action (Wilson or 
staggered) on the finest lattice and incorporate /i by this rule, which we call the standard 
procedure. In coordinate space, it amounts to the substitutions 

"^{x, t) = e^^'^ix, t) , '^^(f , t) = e-'^*^(f , t) . (2) 

Hence /i is treated on the lattice consistently as an imaginary constant Abelian gauge 
potential A4 = ifi. Therefore, and are "parallel transported to t = 0", which 
confirms the gauge invariance of this procedure. This observation holds for all sorts of 
lattice actions, hence the standard procedure (@) can always be applied. 

We preserve this useful property also under the RGT by using the transformation term 

yp^/M,|r^"A'^'M^^f/'^f/]^ ^f^q,'[x',t') = e^"*'^'(x',t'), "^^'(x',t') = e-'^"*'^'(f',t')-(3) 

If t G [0,/3], then t' G [0,f3/n], and since distances are measured after the RGT in coarse 
lattice units, the transportation distance to t = t' = is divided by n. This is compensated 
by substituting n/i in the coarse fields. 

The result of the RGT is 5"["m^'^"m ^'^ f/']^ and it is easy to see that this is identical to 
the following construction: first let /i = and perform the RGT to arrive at S'[^', ^E^', U']. 
Now we include in the blocked action again by the standard procedure (0) (where /x is 
multiplied by n due to the rescaling to the new lattice units). 



3 



Now the RGT can be iterated. In order to arrive at a chemical potential /i in the final 
action, we start with /i/nN on the finest lattice, and perform RGTs of block factor n, 
in analogy to the mass. The limit nN oo yields a perfect action with /i. 

From the analysis of one RGT we conclude that we can switch off /i in the beginning, 
iterate the RGT all the way to the perfect action, and then include /i again by the 
procedure (Q). This leads to a perfect action at chemical potential /i 

The standard procedure to include fi preserves perfectness, it does not cause any lattice 
artifacts and it is therefore perfect itself We emphasize that this result applies to any 
perfect action for the fully interacting quantum theory. 

In practice, such actions are hard to construct (for recent reviews, see p). What is 
more realistic is a classically perfect approximation, where the functional integral in eq. 
(|l|) is simplified to a minimization of the exponent on the right-hand side. A number of 
2d studies suggest that classically perfect actions also suppress the lattice artifacts very 
strongly. Consideration of the transformation term shows that including /i in a classically 
perfect action by the procedure (0) leads again to a classically perfect action with /i. 

An alternative improvement program for lattice actions starts from a standard action 
and tries to eliminate the artifacts order by order in the lattice spacing {Symanzik's 
program). Assume that this is realized to some order at = 0, and we include fi once 
more by the standard procedure. The above properties suggest that we obtain an action 
with fi, which is still free of lattice artifacts to the same order as it was the case at fi = 0. 

3. Scaling behavior of a free fermion gas as an example 

A perfect lattice action for free fermions with mass m at /x = reads 



where a is a free RGT parameter. At m = 0, a = oo the action is chirally symmetric 
but non-local. For any finite a the action is local: its couplings in coordinate space decay 
exponentially. Since the limit nN oo relates the system directly to the continuum, 
there is still chiral symmetry in the observables 0. Furthermore, the action itself also 
keeps a remnant continuous chiral symmetry [|To| . 

For applications, it is important to optimize the locality, so that the truncation of the 
couplings is not too harmful. This is achieved by choosing a = (e"^ — m — l)/m? p]. 

According to section 2, the corresponding perfect action at finite n is characterize by 
A(p,p4 -|- ifj). A successful and yet applicable approximation is a truncation to the cou- 
plings in a unit hypercube by means of periodic boundary conditions. The corresponding 
couplings at /i = for various masses are given in Ref. |1TT|| . 

The RGT used for the above perfect action for Wilson-type fermions is based on the 
usual block average (BA) scheme. If one constructs the analogous perfect action for 
staggered fermions [|l^ , then the truncation to the same number of degrees of freedom as 
before (couplings distance components ±1, ±3 resp. 0, ±2 at m = 0) is not satisfactory. 
The locality - and hence the quality after truncation - can be improved significantly by 
using instead a blocking scheme that we call "partial decimation" (PD) ||13|| . 

^What has been achieved - to a good accuracy - is the 0{a) improvement of QCD with Wilson fermions 
I?). However, going beyond that is not realistic in the foreseeable future. 
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We now focus on m = and temperature T = 0. The pressure and baryon density of 
the free fermion gas are given by 

" f^B = — o continuum, and obtained from 



d^p detA(p,p4) 1 "5 p fv, 1 ff 

log- — ——; — - , ub = 77"^ -< on the lattice. (5j 



TV {2tc)^ detA(p,p4 + i/i) ' 3 9/i 

For comparison, we insert a number of lattice fermion propagators A. For Wilson-type 
fermions, we probe the Wilson fermion with Wilson parameter r = 1, a Symanzik im- 
proved version thereof called D234 with additional couplings on the axes [Q, and the 
truncated perfect "hypercube fermion" (HF). For staggered fermions we insert the stag- 
gered standard action, the Naik fermion JTSf (which is Symanzik improved in the same 
way as the D234 fermion), and again the truncated perfect actions for the BA and the 
PD scheme. Fig. 1 shows the scaling ratios and ub/^j? for this set of fermions. 

Of course, they all reproduce the correct continuum values at yU — 0. As /i increases, 
the standard actions deviate very soon, the Symanzik improved actions do well up to a 
certain /i and then collapse completely, but the truncated perfect actions keep close to 
the continuum value up to remarkably large /x. Their scaling region is extended by one 
order of magnitude compared to the standard actions. 

One could also consider XbI l^^-, where xb = ^'^s is the baryon number susceptibility, 
but in all these cases the qualitative behavior is very similar, and it is also in agreement 
with the scaling of P/T^ at fi = [^,0, hence it is manifestly systematic. 

4. The fermion density 

We now consider the fermion density n^. Q It is well known that the standard lattice 
actions suffer from an upper bound for nf. We choose the normalization such that naive 
fermions saturate at ri/ = 32. This corresponds to an occupation of each site by 2'^ species, 
each with spin up and down. 

For Wilson fermions the doublers are heavier, hence fi must be larger to exceed their 
masses, and the saturation is delayed. However, for r ^ 1 this is unimportant in the limit 
/i oo, so we have again nj^max = 32. The case of r = 1 is special: here the mass of the 
time-like doublers diverges, hence nj^max = 16. 

From Pauli's principle one might be tempted to conclude that = 32 is an absolute 
upper bound for all lattice actions. However, since this bound is a lattice artifact, this 
leads to a puzzle about the behavior of a perfect action. The solution is that the above 
argument is wrong. Considering 

S[^, = E E ^^,t e^(*-^) A-\x -y,t~s) ^y^,s , (6) 

x,t y,s 

we recognize that the exponential growth of the factor containing fi is limited by the 
maximal coupling distance in the temporal direction, which occurs in A~^. (In fact, such 
a saturation limit also exists for bosons.) For the perfect action this distance is infinite, 
hence n/ has no upper bound, in agreement with the absence of artifacts. To turn it 
the other way round, we can conclude that (in infinite volume) no perfect action with 

•^At this point I thank especially M.-P. Lombardo, who suggested to me to consider also nf{^) itself. 
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Figure 1. The scaling of Wilson-type fermions (on top) and staggered fermions (below), 
both in a number of variants discussed in the text. The D234 and the Naik fermion are 
Symanzik improved, and the fermions denoted as HF, BA and PD are truncated perfect. 

couplings only in a finite range can exist, and this is indeed correct. For the hypercube 
fermion we have again [t — s)max = 1 — * nj^rnax = 32, but the saturation is significantly 
delayed. The D234 action has (t — s)max = 2, so one could expect nj ^as = 64. But there 
is a cancelation going on, similar to the Wilson fermion at r = 1, which reduces Uf^^nax 
again to 32. If we keep that structure and vary the couplings, then Uf does rise up to 64. 
But in both cases, the saturation occurs at practically the same value of /i as it happens 
for the r = 1 Wilson fermion. This behavior is shown in Fig. 2, which also illustrates the 
Wilson fermion at various values of r. For instance at r = 2, nf{fi) rises quite steadily to 
32. If r decreases towards 1, then half of the particles in the spectrum turn very heavy. 
The curve reaches a first plateau at 16, and when /i catches up with this heavy species, 
it performs a second jump up to 32. This looks similar if we approach r = 1 from below. 
If r is exactly 1, then the plateau at 16 extends to infinity. 

We also show these curves for various types of staggered fermions. The standard bound 
is 8, but if couplings over temporal distance 3 are involved, then this bound is amplified 
by a factor of 3. Here we also show some Symanzik improved actions with diagonal 
couplings from Ref. ITO. One type called p4 has maximal temporal coupling distance 



"^The only exception for matter fields is the case d = 1, which does, however, not lead to a contradiction. 
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Figure 2. T/ie fermion density for some lattice fermions mentioned in the text. 

2 and saturates at ri/ = 16, the others couple over distance 3 hke the truncated perfect 
fermions. Also here the saturation is delayed, and we mention that especially the p6 
action is doing very well in the scaling tests too. 

5. Generalizations 

A number of generalizations are straightforward, for instance the steps to finite tem- 
perature and finite mass. The latter has been considered in Ref. [^] for naive fermions, 
and the behavior is similar in all other cases: the saturation is delayed as m grows, but 
the upper bound is not affected. 

It is fashionable to use anisotropic lattices in thermodynamic simulations. In particu- 
lar, ^ = ttspatiai I ^temporal > 1 is uscful, becausc it provides a better resolution of a decay 
in the temporal direction. Remarkably, the couplings in a perfect action are indepen- 
dent of C,- This property still holds after truncation by periodic boundary conditions. 
In the previous plots, we would reproduce the same curves, and the horizontal axes now 
means atemporav 1^'- However, the lattice units are defined by aspatiai, so /x is multiplied by ^. 



6. Conclusions 

So far, /i has only been included in standard lattice actions. We have seen how to do 
this in general. Of interest are applications to improved actions, and there the improved 
quality (perfectness, classical perfectness, or Symanzik improvement to some order) is 
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preserved under the inclusion of /i. 

We have shown that an improved action allows us to go beyond the standard limits for 
the baryon resp. fermion density (in spite of a widespread believe that this is impossible). 

Applicable truncated perfect actions extend the scaling region of a free fermion gas by 
one order of magnitude, which suggests that simulations with them can be carried out 
on very coarse lattices. This is not a direct remedy of the sign problem mentioned in 
the introduction (because that is a problem of statistics, whereas the improved action 
cures systematic errors), but it helps indirectly, because larger physical volumes can be 
handled. 

For the fermion-gauge couplings and for the pure gauge part, there is much work going 
on to construct Symanzik improved or approximately perfect actions If this is 

achieved at n = 0, then the extension to /i 7^ is solved by the prescription shown here. 

Hopefully this is a step towards more conclusive simulations for ub ^ 0. 

Most of this talk is based on Ref. j^; I am indebted to U.-J. Wiese for his collaboration. 
I also thank M.-P. Lombardo and F. Karsch, both, for organizing this workshop and for 
instructive discussions. 
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